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Abstract 

We analyze the critical behavior of the dephasing rate induced by short-range 
electron-electron interaction near an Anderson transition of metal-insulator or quan- 
tum Hall type. The corresponding exponent characterizes the scaling of the transi- 
tion width with temperature. Assuming no spin degeneracy, the critical behavior can 
be studied by performing the scaling analysis in the vicinity of the non-interacting 
fixed point, since the latter is stable with respect to the interaction. We combine 
an analytical treatment (that includes the identification of operators responsible for 
dephasing in the formalism of the non-linear sigma-model and the corresponding 
renormalization-group analysis in 2 + e dimensions) with numerical simulations on 
the Chalker-Coddington network model of the quantum Hall transition. Finally, we 
discuss the current understanding of the Coulomb interaction case and the available 
experimental data. 
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1 Introduction 



Localization-delocalization quantum phase transitions form a broad and ac- 
tively developing field of condensed matter physics, see Ref. [1] for a recent 
review. In particular, the electron-electron interaction effects at the transition 
represent one of central research directions [2,3,4,5,6,7]. 

Quite generally, the impact of interaction onto low-temperature transport and 
localization in disordered electronic systems can be subdivided into effects of 
(i) renormalization and (ii) dephasing. The renormalization effects, whose role 
in diffusive systems was investigated by Altshuler and Aronov [2], are governed 
by virtual processes and become increasingly more pronounced with lowering 
temperature. Finkelstein developed a renormalization-group (RG) approach 
[3] in the framework of the non-linear a-model in order to treat these effects 
together with localization phenomena. More recently, this research direction 
attracted a great deal of attention in connection with experiments on high- 
mobility low-density electronic structures (Si MOSFETs) giving an evidence 
in favor of a metal-insulator transition [6,7,8]. It was shown that this transi- 
tion can be explained in the framework of the cr-model RG for a system with 
N > 1 valleys [9] (formally N should be large but in practice N = 2, as in 
Si, is already sufficient). Indeed, a detailed analysis has confirmed that the 
RG theory describes well the experimental data up to lowest accessible tem- 
peratures [10,11]. Very recently, it was shown [12] that interaction effects are 
also of crucial importance in topological insulators, where they induce novel 
critical states. Another long-standing issue in the field concerns the interplay 
of interaction and multifractality of wave functions at Anderson transitions. 
This problem has been recently addressed in the context of superconductivity 
near the mobility edge. Focussing on the short-range interaction in the Cooper 
channel, Ref. [13] predicted that the multifractality of critical wave functions 
strongly enhances the superconducting pairing correlations. 

Dephasing effects are governed by inelastic processes of electron-electron scat- 
tering at finite temperature T. The dephasing has been studied in great detail 
for metallic systems where it provides a cutoff for weak-localization effects 
[2]. As to Anderson localization transitions, dephasing leads to their smear- 
ing at finite temperature. The dephasing-induced transition width scales as a 
power-law function of T. In the case of quantum Hall transition, this scaling 
of the transition width, AB oc T K , has been experimentally explored in many 
works. While the value of k was a matter of controversy, there seems to be 
a consensus now that the results k = 0.42 ± 0.04 [14], k = 0.42 ± 0.01 [15] 
properly characterize the quantum Hall transition (when it is not masked by 
macroscopic inhomogeneities) . Scaling near the transition with varying tem- 
perature was also experimentally studied at the 3D Anderson transition in 
doped semiconductors [16,17]. 
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In this paper we study the scaling properties of dephasing rate at criticality. 
We consider the situation when the time-reversal invariance is broken, so that 
the system belongs to the so-called unitary symmetry class. In particular, the 
quantum Hall transition belongs to this symmetry class. We further assume, 
following Refs. [18,19] that the interaction is of short-range character. This 
greatly simplifies the analysis, since in this case the interaction is irrelevant 
in the RG sense [4,18,20], so that all scaling properties at the transition are 
governed by the RG flow near the non- interacting fixed point. The goals of 
this work follows. First, we identify the operators controlling the de- 

phasing within the framework of the non-linear cx-model. Second, we perform 
the RG analysis of these operators for the cr-model in 2 + e dimensions. Third, 
we present a numerical investigation of the corresponding wave function cor- 
relation functions by using the Chalker-Coddington network as a a model of 
the quantum Hall transition. Finally, we make conclusions concerning the ex- 
ponent k governing the temperature scaling of the transition and discuss the 
current understanding of the Coulomb interaction case, as well as relations 
between the theoretical and experimental findings. 



2 Dephasing rate and exact eigenfunction of the non-interacting 
problem 



2.1 Model 



We consider a disordered interacting electronic system characterized by the 
Hamiltonian 

H = Ho + Hint ■ (1) 

The one-particle term in the Hamiltonian, 



H = J drm\r) 



2m P 



-iV - A(r)) 2 + V dis (r) 



tt(r), (2) 



is assumed to describe a system with broken time-reversal invariance (unitary 
symmetry class) at criticality. This may be a 2D system near quantum Hall 
transition (in which case B = rot A is a uniform magnetic field), or a system in 
d > 2 dimensions at the Anderson transition point (with a uniform or random 
magnetic field B(r) breaking the time reversal invariance). The interaction 
part of the Hamiltonian % has the form 

%nt =\J d ri dr 2 U{r x - r 2 )ttt( ri )*t( r2 )tf ( r2 )* ( ri ) , (3) 

with a short-range interaction potential U(r). 
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Fig. 1. The first-order (Hartree and Fock) interaction corrections to the thermo- 
dynamic potential. The wavy line denotes the interaction and the solid line the 
electron Green's function. 

2.2 First-order correction to the thermodynamic potential 



Following Ref. [18], we begin by considering the first order interaction correc- 
tion to the thermodynamic potential. The corresponding diagrams are shown 
in Fig. 1; their total contribution reads 

= ^E^^KM / dr 1 dr 2 U(r 1 -r 2 ) B a p{r 1 ,r 2 ) , (4) 

aft 

where n/(e) is the Fermi-Dirac distribution, <f) a ( r ) and e Q are exact eigenfunc- 
tions and eigenenergies of the non- interacting Hamiltonian (2), and 

# Q/3 (ri, r 2 ) = <p a {ri)<pp{r2) - a (r 2 )^(ri). (5) 

After averaging over disorder, the contribution (4) becomes 

{Q W) = J ^^ nf (E)n f (E + u) J drtdrzUin-r^Jdir^^E^), (6) 

where A = 1/ Vd,L d stands for the mean-level spacing, L is the system size, and 
Vd the density of states. The function 

£i = ^" E (l-Mrx, r 2 )\ 2 5(E + u- e a )5(E - e p )\ (7) 

1 a/3 ^ ' 

describes correlations of two eigenstates of the non-interacting Hamilto- 
nian (2). In general, for the system size L — > oo, the function /Ci exhibits 
the following scaling behavior 

id(\ r i ~ r 2\Y 2 s. (\ri- r 2 \\ 



x < 1, 

X > 1 
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where L w = L(A/|w|) 1//d is the length scale set by the frequency u. For two 
adjacent in energy eigenstates one has u ~ A and L u ~ L, for larger frequen- 
cies the scale is much less than L. It is important that the exponent is 
positive, fi2 > 0, due to Hartree-Fock antisymmetrization of wave functions in 
Eq. (5). This means that for distances shorter than L u the correlation function 
(8) is suppressed at criticality as compared to the metallic system. This should 
be contrasted to multifractal correlations of wave functions (without antisym- 
metrization) that are enhanced at criticality. We will discuss this point and, 
more generally, the spectrum of critical exponents within the a-model frame- 
work in Sec. 3. 

We assume the electron-electron interaction of the following form: 



U(R) = u 



1 + (R/a) 



-A/2 



(9) 



with A > d. This yields 



<n«> 



dEdu 
A 



n f (E)n f (E + u)u(L u 



(10) 



where 



u{L u 



d + /i 2 < A, 
A = d + nz, 
d < A < d + /i 2 . 



Following Ref. [18], we can consider u(L w ) as renormalized interaction param- 
eter. It is seen that for A > d the electron-electron interaction is indeed irrele- 
vant near non-interacting fixed point u — 0. While the above calculations are 
based on evaluation of the Hartree-Fock contribution to the thermodynamic 
potential, this conclusion is in agreement with the analysis of Ref. [20] based 
on Finkelstein non-linear cr-model. 



2.3 Dephasing rate: Second-order correction to the electron self-energy 



In order to compute the dephasing rate (or, more accurately, the out-scattering 
rate which in this case coincides with the dephasing rate), one needs to evaluate 
the second-order interaction correction to the imaginary part of the self-energy 
(see Fig. 2). In spirit of Ref. [21], we define the self-energy of a given single- 
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Fig. 2. Second-order interaction contributions to the electron self-energy that de- 
termine the dephasing rate. 

particle state a, 

Sa(e) = g J dr 1 dr 2 dr 3 dr 4 U(r 1 - r 2 )U(r 3 - r 4 )B* ol p(r x ,r 2 )B^{r 1 ,r 2 ) 

n f (ep)[l -ri/(e 7 )] + ra / (e*) [n f (e 7 ) - n f {e p )\ 



xB* 5 (r 3 ,r 4 )Bp a (r 3 ,r 4 )J2 

^5 e + e/3 - e 7 - e 5 + zO 



(12) 



The next step is the evaluation of the imaginary part of the averaged self- 
energy 

Z R (E,e) = A(^Z*(e)5(E-e a )Y (13) 

Since we consider the situation at criticality, where fluctuations may be strong, 
the averaging is not always an innocent procedure. We will return to the 
justification of the averaging in Sec. 4.4. 



The result of averaging can be presented in the form 
lmJ: R (E,e) = -it 



x [1 - n f (e')] + [n f (e') - n f (e' + Q)]n f (e + fi)|/C 2 ({r i }, E, e, e\ ft) , 

(14) 

where the correlation function K, 2 ({rj}, E, e, e', CI) is defined as follows 

A 4 / 

fc2({rj},E,e,e',n) = —( ^ B* ap {r ll r2)B h {r ll r 2 )B* l5 {r 3 ,r 4 )Bp a {r 3 ,r 4 ) 
x6(E - e a )5{e' + n- ep)8{e' - e 7 )«5(e + Q- e s )). (15) 



To determine the scaling of the dephasing rate, we now specify characteristic 
values of energy variables. First of all, we are interested in dephasing at the 
mass shell (E = e) and at characteristic energy E ~ T '. Since we are not 
aiming at calculating the numerical prefactor of order unity, we can simply 
set E = 0. Second, we will see that the characteristic values of the integral 
variables e' and Q are set by the temperature. Performing integration of the 
Fermi- function factor over energy e', we find 



ImS K (0,0) 



"2^3 (A/ dr ^j U ^ - r ^ U (rs ~ r 4 ) / 



dQQ 



x 



{coth ^ - tanh ^JC 2 ({ rj }, 0, 0, e' ~ T, Q) . 



(16) 



In order to proceed further, we need to know the scaling behavior of the 
function /C 2 . Assuming that e' ~ (this is sufficient for our purposes, as both 
these frequencies are of order of temperature), we have for |r x — r 2 |, |r 3 — r 4 | <C 
R ^ Lq the scaling form [18] 

where R = (r\ + r 2 — r% — T4)/2. The scaling behavior (17) can be motivated 
as follows. At distances R ~ the correlations between wave functions 
at r 12 , on one hand, and r 34 , on the other hand, decouple. Thus, scaling 
of /C 2 reduces to that of a product of two independent correlators /Ci, i.e., 
(|t*i — T 2 |/-R) M2 (|r 3 — r 4 |)/i?) M2 . The exponent a describes the scaling with 
respect to a remaining scaling variable R/Lq. For R ^> Lq the correlations 
quickly (exponentially) decay, so that this range of R is not important for the 
integral. 

Combining Eqs. (16) and (17), we obtain 

|ImS fl (0,0)| ~ u d T [ T dil dRu 2 (R) (^-) . (18) 

JO J a \LqJ 

The result of integration over R depends on the relations between A, /i 2 , a 
and d. Let us first assume that a > —d. Then, if A > d + /i 2 , we find 

( {T/5) 1+2 ^/ d , d - 2/i 2 + a > 0, 

— oc 5u 2 (a) | {T/5) 1+2 ^l d \n(5/T), d - 2/i 2 + a = 0, (19) 
Tlf [(T/6) 2+a / d , rf-2/i 2 + «<0. 
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In the case A = d + /i 2 , we obtain from Eq. (18) 

C {T/5) 1+2 ^/ d ln 2 (5/T), d - 2fi 2 + a > 0, 

— oc 5u 2 (a) | (T/5) 1+ W d ln 3 (<VT), d - 2^ 2 + a = 0, (20) 
Ttp [{T/5) 2+a / d , d-2/i 2 + a<0. 

Finally, if A < d + /i 2 , we get 

f(r/5)- 1+2A ^, 3d-2A + a>0, 

— ot5u 2 (a) I (T/5y 1+2X / d \n(5/T) 1 3d - 2A + a = 0, (21) 
^ [(T/5) 2+Q / d , 3d-2A + a<0. 

Here 5=1/ i>da d is the ultraviolet energy cutoff. 

For a < —d, the integration over Q is dominated by the infrared cutoff which 
should be chosen self-consistently at 1/t v , yielding 

J_ \Tu 2 (a)\ In u(a) |, a = -d, 

K 1 5 [Tn 2 (a)/5]^ d/Q , a < -tZ. 1 J 

Thus, the temperature behavior of the dephasing rate depends on the expo- 
nents fi2 and a as well as on dimensionality d and the index A characterizing 
the decay of interaction. The exponents /i 2 and a belong to a set of exponents 
characterizing correlations of wave functions at criticality; their values depend 
on the particular critical point under consideration. In the next section we rep- 
resent the correlation function /C 2 in terms of the operators in the non-linear 
sigma model. This will allow us to substantiate the above scaling arguments 
by the RG analysis and to evaluate the necessary exponents within the 2 + e 
expansion. 



3 Field theoretical approach to wave function correlations 

3. 1 Non-linear sigma model 

The effective theory characterizing long-distance low-energy physics of a dis- 
ordered electronic system is the non- linear cr-model [22]. We will use its replica 
version; the same calculations can be performed in the supersymmetric formu- 
lation. The field variable of the theory is the matrix field Q(r) of size 2n x 2n 
that obeys the non-linear constraint Q 2 (r) = 1. The non-linear sigma model 
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action has the form 



S a [Q] = - 9 -Jdrti(VQ) 2 + 9 -^- J drtiQA, (23) 

where 

A= I 




The parameter g in front of the kinetic term is identified with the dimension- 
less longitudinal conductance (in units of e 2 /h). In the RG framework, it is 
interpreted as a running coupling constant of the theory. The second term 
in Eq. (23), which breaks the U(2n) symmetry down to U(n) x U(n), nor- 
mally has an imaginary prefactor proportional to frequency e. We choose the 
constant in front of this term to be real, which is more convenient for RG 
analysis. 

In the case of a 2D system in transverse magnetic field, there is an additional 
contribution to the action that has a form of the topological term [23] 

Stop [Q] = / dr tr QW x QW y Q . (25) 

The angle 6 in front of the topological term is given by the fractional part of 
the 2ng xy , where g xy stands for the dimensionless Hall conductivity (in units 
e 2 /h) [24]. 

We now proceed by expressing the correlation function /C 2 in terms of eigen- 
operators of the cr-model. This calculation is only based on the symmetry of 
the theory, so it is equally applicable to the quantum Hall transition in 2D 
and to the Anderson transition in d > 2 dimensions. 



3.2 Expression for K.2 in terms of the eigenoperators 



It is convenient to express first the function /C 2 in terms of the exact single- 
particle Green's function G^' A (r,r f ): 



7T J 



JC 2 ({ rj },E,e,e',n) = — (ImG^rOlmG^nfo,^ 



lmG R 



x ImGp(r 2 , r 3 ) - ImG^ +n (ri, r 3 ) ImG^(r 2 , r 4 



Ri 



£ +n{r u r 4 ) 
(26) 



In order to extract the exponent a, it is enough to consider scaling of 
K,2({rj}, 0, 0, 0, 0) with the system size L for fixed \rj — r^j. Then, following 
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standard steps (see e.g. Ref. [25]), we obtain 



JC 2 ({ rj }, 0,0, 0,0) 



(2tt7) 4 



tr 



AQafe(r-i)AQ 6 a(r 4 ; 



tr 



AQcd(r 2 )AQ dc {r 3 ) 



+ tr 



AQab(ri)AQU r 4) A Qcd(r 2 )AQ da (r 3 ] 



(27) 



Here 7 = (irud)" 1 , the symbol tr stands for the trace over retarded-advanced 
space only and (...) denotes the averaging with the a- model action S a . It is 
important to stress that the replica indices a, b, c, d are fixed and different: a 7^ 
b 7^ c 7^ d. It reflects the fact that we are dealing with four different eigenstates. 
When deriving Eq. (27), we have assumed that all distances \rj — exceed 
the mean free path, so that the averaged Green's functions between any two 
points can be neglected. (This assumption does not affect the scaling but 
simplifies calculations.) On the other hand, we can think about all points as 
located close to each other from the point of view of the a-model (say, on 
the scale of several mean free paths), so that we omit space argument of Q 
matrices in what follows. 



In view of the presence of the fixed replica indices, Eq. (27) is not convenient 
for further computations with the action S a which is U(n) x U(n) invariant. 
In order to obtain a more convenient, U(n) x U(n) invariant expression, we 
use the fact that the correlation function is invariant with respect to global 
unitary rotations in the replica space. This allows us to average the expression 
in angular brackets in Eq. (27) over such global unitary rotations. In general, 
the U(n) x U (n) invariant operators of the kth order in Q matrices have the 
form 0\ = Tr(AQ) fcl . . . Tr(AQ) fcm , where A = {k\, . . . , k m } is a partition of 
the number k = k\ + . . . + k m , such that k± > k 2 . . . > k m . In particular, for 
k = 1 we have one operator, {1}, for k = 2 two operators {2} and {1, 1}, for 
k = 3 three operators {3}, {2, 1}, and {1,1,1}, for k = 4 five operators {4}, 
{3, 1}, {2, 2}, {2, 1, 1}, and {1, 1, 1, 1}, and so on. 



Since the correlator /C2, Eq. (27), contains four Q matrices, after the averaging 
over U(n) x U(n) global rotations the result is expressed in terms of the 
invariant operators of even orders < 4, i.e. with k = 4 and k = 2 (and in 
addition a constant term corresponding to k = 0). Specifically, we obtain (see 
Appendix B) 

K = j0ff{g c MQ\ + G\ , (28) 

where operators Oj[Q] and corresponding coefficients are given as 
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4 [Q]=Tr(AQ) 4 , (29) 

3 ,i[Q]=Tr(AQ) 3 TrAQ, (30) 

2 , 2 [g]=Tr(Ag) 2 Tr(Ag) 2 , (31) 

2 , 1 , 1 [g]=Tr(Ag) 2 (TrAg) 2 , (32) 

Oi,i,i,i[g] = (TrAg) 4 , (33) 

2 [g]=Tr(Ag) 2 , (34) 

L1 [g] = (TrAg) 2 (35) 



and 



a 
c. 

c 

^2,1,1 



3,1 — 



2,2 — 



c x 



-3 + 13n + 16n 2 + 4n 3 
4n 2 (-l + n)(l + n) 2 (2 + n)(3 + n) 
-7 + 3n + 12n 2 + 4n 3 



n 2 ) 2 (3 + n)(2 + n) 

,2 i 



2n 2 (-l - 
-3 - 21n + 20n 2 + 32n 3 + 8n 4 
> 2 (-l +n 2 ) 2 (3 + n)(2 + n) ' 
3 + 2n 



2n 2 



■n)(l + n) 2 (3 + n) 



1,1,1 — 



C 2 = 



CVi = 



x -1 

5 + 5n + 2n 2 
8n 2 (-l + n 2 ) 2 (3 + n)(2 + n)' 

-3 + 3n + 4n 2 
2n 2 (-l + n) 2 (3 + n) ' 
3 + 5n + 4n 2 



2ra 2 (-l + n 2 ) 2 (3 + n)(2 + n)' 
76n 3 - 24n 4 



-9 + 93n + 4n 2 
°~ 2ra(-l +n 2 ) 2 (3 + n)(2 + n 



(36) 
(37) 
(38) 
(39) 
(40) 
(41) 
(42) 
(43) 



While the operators Oj[Q] are explicitly invariant with respect to the U(n) x 
U(n) symmetry group of the cr-model, they are not eigenoperators of the renor- 
malization group. The eigenoperators that are linear combinations of Oj[Q] are 
in fact fully determined by the symmetry group (and realize its different repre- 
sentations, similarly to conventional spheric function realizing representations 
of the rotation group) [29]. The seven RG eigenoperators can be enumerated 
by Young frames: P 4 , P 3 i, P 2j2 , P2,i,i, Pi,i,i,i, P2, and P^i [30,31]. To express 
these RG eigenoperators in terms of basis operators Oj[g] one can use the 
results of Refs. [31,32] or the one-loop renormalization [30,26]. Finally, using 
Eq. (28) we represent the correlation function /C 2 as a linear combination of 
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the eigenoperators of RG, 



A 4 / 4n 2 4n 2 + 8n + 3 

2 (2tt 7 ) 4 \4n 2 (n 2 - l)2^ 2 ' 2 ^ 6 n 2 (n + l) 2 (n 2 + n - 2) 2 ' 1>llVJ 

4n 2 + 16n + 15 n r _..\ fAt . 
+ 2 n »(B + l)»(n + 2)(n + 3) ^-" W1 /- (44) 

It is remarkable that only three eigenoperators (out of seven) enter the ob- 
tained expression for the correlation function /C 2 . The eigenoperators P4, P^x 
and P 2 with leading singularities in the replica limit n = do not contribute 
to /C 2 . 



4 Temperature dependence of dephasing rate at criticality 



4-1 Anderson transition in 2 + e dimensions 



In d = 2 + e dimensions with small e the Anderson transition takes place 
at large dimensionless conductance g. This allows one to explore the critical 
behavior within the e-expansion. The /3-function governing the renormalization 
of the conductance is known up to the five-loop order [27,28] 

dt 

P{t) =et- 2t 3 - 6t 5 + 0(t 7 ) , (45) 



dlnL 



where t = l/2irg. The metal-insulator transition occurs at the critical point 
defined by /3(t*) = that yields 



1/2 3 / e x3/2 



.2/ 2 \2J 
The localization length index v reads 



-) +0(e 5 / 2 ). (46) 



1/ = -!//?'(**) = 1 - ^ + 0(6) . (47) 



Further, according to Ref. [31], the anomalous dimensions of the eigenopera- 
tors Pj are given up to four-loop order as 

J Pj (t) = a 2 p(t) + C(3) C3 t 4 + 0(t 5 ), pit) = t+-(n 2 + l)t 3 + -(n 2 + 7)t 4 (48) 

where coefficients a 2 and C3 are summarized in Table 1. It is worthwhile to 
mention that the coefficient C3 is proportional to a 2 . Therefore, in the replica 
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Table 1 

Coefficients of 7-functions from Ref. [31]. 





«2 


C3 


Pa 


12 - 8n 


12(3-2n)(2-n)(3-n) 


p 3,i 


4 - 8n 


4(1 - 2n)(6 - 13n + n 2 ) 


p2,2 


-8n 


84re 2 


^2,1,1 


-4 - 8n 


4(1 + 2n)(6 + 13n + n 2 ) 


-Pi, 1,1,1 


-12 -8n 


12(3 + 2n)(2 + n)(3 + n) 


P 2 


2 - 4n 


2(l-2n)(l-n)(3-n) 




-2 - 4n 


2(l + 2n)(l + n)(3 + n) 



limit n — > 0, the anomalous dimensions of the eigenoperators involved in 
Eq. (44) become 

7P 2 , 2 (t) = 0(t 5 ), (49) 
7jv,i(*) = - 4t - 6t3 + 24C(3)t 4 + 0(t 5 ), (50) 
7Pl.i,i.i (*) = ~ m - 28t3 + 216C(3)t 4 + 0(f). (51) 

The exponent \ii is determined by the anomalous dimension of the eigenoper- 
ator Pi 1 [30,26], see Appendix A. In the replica limit n — > 0, 

7Pia (t) = -2t - 3t 3 + 6C(3)t 4 + 0(t 5 ). (52) 

Using Eq. (46), we find the exponent \ii in 2 + e dimensions [26,31], 

M2 = -7p m (^) = v/27- ^C(3)6 2 + 0(e 5 / 2 ). (53) 

The exponent a is determined by the maximal value of the anomalous dimen- 
sions (with oppposite sign) of operators Pi,i, ^2,1,1 and -Pi, 1,1,1 at the critical 
point: 

a = max{-7p 22 (t*), -7^,1,1 (**), -7ft,i,i,i (**)}■ ( 54 ) 
The anomalous dimensions of the operators P*2,i,i and -Pi, 1,1,1 are negative, 
whereas 7p 2 2 vanishes within the known accuracy: 

a = 0(e 5 / 2 ). (55) 

Therefore, in the case A > d + y2e (for brevity we keep here only the leading, 
one-loop, contribution to /i 2 ), the temperature behavior of the dephasing rate 
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is given by 



1 /T\ 

-oc 5u 2 {a) [=-) . (56) 

T<f> V J 



In the special case A = d + v^e, 



— oc 5u 2 (a) (^-] +VT \n 2 {5/T). (57) 



Only in the regime d < A < <i+v 2e, the temperature behavior of the dephasing 
rate is determined by the exponent A: 

— oc5u 2 (a)(^) . (58) 



Let us focus on the case of the "most short-range" interaction, A > d + /i2, 
when t^ 1 oc T p with p = 1 + 2fi 2 /d. The scaling of the dephasing length is 
then given by 

The exponent belongs to a class of dynamical critical exponents, as it 
governs the scaling of a characteristic length scale (dephasing length) with 
a variable having the dimension of energy. It is worth emphasizing that the 
localization probems possess rich physics and are in general characterized by 
several dynamical exponents. This fact has not always been appreciated in 
the literature. In the present case, one should distinguish the exponent z T 
controlling the scaling with temperature from the exponent z governing the 
scaling with frequency. The latter exponent has a trivial (non-interacting) 
value z — d for the short-range interaction. 

The transition width induced by inelastic scattering scales as T K , where the 
exponent k is found by comparison of the dephasing length (59) with the 
localization (correlation) length £ oc \P — P c \~ u , where P is the parameter, 
driving the transition (e.g., electron concentration or disorder strength). This 
yields 

i i + Wj 

k= = . (60) 

ztv vd 

Substituting the above formulas for fi 2 and u, we get the following results for 
the dynamical exponent zt and for the index k up to four-loop order of the 
e-expansion: 



z T = 2- 2v / 2e 1 / 2 + 5e - 4 v / 2e 3/2 + 0(e 2 ) ; (61) 
K = e + v / 2e 3/2 + e 2 + e 5/2 /v / 2 + 0(e 3 ). (62) 
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Fig. 3. Hartree and Fock contributions to the wave function correlator JCi as denned 
in Eq. (7) for pairs of wavefunctions neighboring in energy. It is seen that in the 
scaling regime both correlation functions follow the same power law, including the 
exponent and the amplitude. The corresponding exponent A2 — —0.52 is known 
from the multifractal analysis [1]. 



4-2 Integer quantum Hall transition 



We have employed numerical techniques in order to calculate K,\^ and the 
corresponding exponents for the quantum Hall transition. Our computations 
are based on the Chalker-Coddington network. We extract eigenvalues and 
eigenfunctions near zero (pseudo) energy from large square systems. In the 
present context the system size L is parametrized with the number of links 
N = 512, 768, 1024 in each direction. Typically of the order of 10 4 wavefunc- 
tions enter the ensemble average for the correlators /Ci^. 

In Fig. 3 we show separately two correlation functions (Hartree and Fock 
terms), the difference of which constitutes the wave function correlator /Ci 
defined in Eq. (7). It is nicely seen that in the scaling regime of point separation 
|ri — r 2 | much smaller than the system size N both correlation functions 
follow the same power law. The corresponding exponent is well known from 
the multifractal analysis of moments of wave functions, A2 — —0.52. It is 
important that not only the exponent but also the prefactor is the same. This 
ensures that the difference between the two correlation functions scales with 
another (subleading) exponent /i 2 . 

Figure 4 shows the correlator /C x representing the difference of the two func- 
tions shown in the previous plot, Fig. 3. As emphasized above, the leading 
power-law contributions to the Hartree and Fock terms cancel, so that /Q is 
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Fig. 4. Correlator K\ representing the difference of the two functions shown in Fig. 
3. The leading power laws cancel and K,\ is determined by the subleading contribu- 
tions. Two types of deviations from pure power-law behavior (straight line in the 
double-log scale) are seen. At small distances the data collapse is not perfect due to 
corrections to scaling originating from the ultraviolet cutoff scale (lattice constant). 
At large distances deviations from the power-law scaling are caused by |ri — T2\ 
approaching N. The extracted value of the power-law exponent is pi — 0.62 ± 0.05. 



determined by subleading contributions. For a pure power-law scaling behavior 
we would have a straight line in the double logarithmic scale. We observe two 
types of deviations from this behavior. At large distances (\ri — r 2 \/N > 0.1) 
there is a considerable curvature which is related to higher subleading terms. 
At small distances the data collapse is not perfect due to deviations from 
scaling related to the ultraviolet cutoff scale (lattice constant) a which leads 
to emergence of an additional scaling parameter a/\r-i — Y2,\. The numerical 
analysis yields the power-law exponent ji 2 — 0.62 ± 0.05. 



In Fig. 5 we show the correlation function K. 2 at fixed small distance between 
pairs of the points, |r x — r 2 | = |r 3 — r 4 | = p, as a function of the distance 
R between the pairs. When the system size N increases, the data approach a 
straight line, corresponding to a power-law dependence on R with the exponent 
~ 1.25. This exponent is equal to 2/i2 within the uncertainty of our numerical 
analysis. Therefore, within our accuracy the exponent a denned in Eq. (17) is 
indistinguishable from zero. 

The above result on the exponent a is further supported by Fig. 6 where 
the correlation function /C 2 is plotted as a function of R/N for fixed p/R. 
According to Eq. (17), this is a direct way to determine the exponent a. We 
see that at R/N < 0.1 the plots are almost flat, which implies a ~ 0. At small 
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0.0001 



Fig. 5. Correlator JC2 at fixed small distances \r\ — r2 1 = |ra — as a function of the 
distance R between the pairs r\ s 2 and ^4. Evolution of /C2 with system size N is 
shown for a quartett of neighboring energies. The average is over 10 6 samples with 
two quartetts selected for each sample; the error bars give one standard deviation 
to indicate the residual statistical uncertainty. It is seen that with increasing N 
the data approach the straight line, corresponding to a power-law dependence. The 
corresponding exponent is ~ 1.25, i.e. equal to 2p2 within the numerical uncertainty, 
implying that a — 0. 

R we observe again the deviations from scaling controlled by the parameter 
a/p, as demonstrated in the inset. 

To make the scaling properties of JC2 particularly clear, we replot in the left 
panel of Fig. 7 the data of Fig. 6, multiplying each trace by (p/R)~ 2tJ-2 with 
2/i2 = 1.25. It is seen that all data (full symbols) collapse on a single-parameter 
scaling curve showing the dependence of K.2{p/R)~ 2 ^ L2 on R/N. This confirms 
that the value of /12 is correct. The open symbols show the data points for three 
smallest values of p\ they deviate from the scaling curve due to corrections in 
a/p. The plateau of the single-parameter scaling curve at small values of R/N 
yields a ~ —0.05 ±0.1 (see the inset in the left panel of Fig. 7 where the same 
data are shown on a log- log scale). The single-parameter scaling curve is 
also shown in the right panel of Fig. 7 on the double-linear scale. 

To summarize, we have found the values P2 — 0.62 ±0.05 and a ~ —0.05 ±0.1 



1 Extracting the bounds on possible values for a from our numerical data is com- 
plicated by the fact that, within the system sizes available to us, corrections in a/p 
and R/L are not negligible, see inset of Fig. 6 and the discussion at the end of this 
section 4.2. The uncertainty in numerical determination of P2 enters the scaling 
analysis (Fig. 7) as an additional source of uncertainty in a. 
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Fig. 6. Scaling behavior of IC2 as a function of R/N for a fixed distance ratio, 
p/R = 1/2, 1/4, 1/8, 1/16, 1/32. Different colors correspond to different system sizes 
(black: 512, red: 768, blue: 1024). According to Eq. (17), we expect the scaling 
cx (R/N) a for R/N <C 1. This scaling behavior is represented by a plateau at inter- 
mediate values of R/N, yielding a close to zero. At larger R/N we see deviations 
from a straight line, since other contributions (scaling with higher irrelevant expo- 
nents) become significant. At small R deviations from the data collapse are due to 
corrections in a/ p. Inset: same data replotted as a function of R/a. At small dis- 
tances no A-dependence is observed. This confirms that deviations from a simple 
power-law scaling (straight line) are controlled by a/ p. 



for the critical exponents describing the correlation functions of interest. These 
values nicely agree with those previously obtained by Lee and Wang [18]. It 
should be stressed, however, that we used systems of a linear size an order of 
magnitude larger than in Ref . [18] . The large system sizes N > 500 were cru- 
cially important for our scaling analysis shown in Fig. 7. Indeed, it is seen there 
that the scaling function reaches the expected power-law behavior oc (R/N) a 
(with a close to zero) only at relatively small values of the scaling argument, 
R/N < 0.05. On the other hand, the scale p cannot be too small, since the cor- 
rections controlled by the parameter a/p become substantial unless p/a > 4. 
Finally, we have a condition R^> p, which in practice requires that R/p > 2. 
Combining all this, we see that in order to have a window of good power-law 
scaling, we need considerably larger than 2 x 4 x 20 = 160. Our system 
sizes reasonably satisfy this requirement. As a result, we get a window of R/N 
approximately between 0.007 and 0.05, where the required scaling takes place. 
On the other hand, relatively small system sizes did not allow the authors of 
Ref. [18] to obtain this scaling window. Indeed, their Fig. 2 shows a function 
that changes by several orders of magnitude without any developed saturation 
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Fig. 7. Left panel: Scaling of (p / R)~ 2f12 K,2 assuming a critical index p2 = 0.625. 
Different symbols correspond to different values of the ratio p/R; different colors 
to different system sizes N, as in Fig. 6. Deviations from the single parameter scal- 
ing for smallest p (open symbols) are due to corrections controlled by a/ p. Inset: 
Same data on a log-log scale. Solid line (corresponding to a = 0) and dotted line 
(corresponding to a = —0.1) indicate the power-law scaling region and show the 
uncertainty range for the numerical value of the exponent a. Right panel: Single- 
parameter scaling function (as obtained from the data points shown by full symbols 
in the left panel) on the double-linear scale. 

plateau. While the authors of Ref. [18] argued correctly in favor of small a, it 
seems difficult to make a reliable conclusion concerning a with only data for 
small systems (as shown in their Fig. 2) at hand. 

Using the obtained values of [i? and a, we can calculate the exponents con- 
trolling the interaction effects, see Eqs. (59) and (60). Assuming the case of 
"most short-range interaction", A > 2 + /i 2 — 2.62, we get 



p ~ 1.62 ; z T ~ 1.23 



K 



0.346 . 



(63) 



While calculating k in Eq. (63), we used the value of the localization length 
exponent, v ~ 2.35 found by Huckestein and coauthors (see the review [33]) 
and confirmed by several later works. Recently, Slevin and Ohtsuki [34] recon- 
sidered the problem and concluded that corrections to scaling are much larger 
that was previously thought. As a result, they obtained a larger value of the 
localization length critical exponent, v ~ 2.59. If this value as used, we get a 
somewhat smaller result for the exponent k, 



k ~ 0.314 



(64) 



We will compare these values to existing experimental results in Sec. 5, where 
we will also discuss expected modifications in the case of Coulomb interaction. 
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4-3 Anderson transition in d = 3 and higher dimensions 

For the Anderson transition in d = 3 the e-expansion cannot give quantita- 
tively reliable predictions for critical exponents. Therefore, one has to rely 
on numerical simulations. The critical exponent of the localization length for 
the unitary symmetry class was found to be v = 1.43 ± 0.04 [35]. While the 
most relevant multifractal exponents (characterizing the scaling of moments 
of wave function amplitudes) A q have been extensively studied, no numerical 
analysis of the subleading exponents, in partiuclar /i 2 and a, has been done to 
the best of our knowledge. This is an interesting direction for future research. 
In particular, an intriguing question is whether the condition d — 2^2 + ot > 
becomes violated in 3D (or, more generally at sufficiently high dimensional- 
ity). According to Eq. (19), this would result in a change of the behavior of 
the dephasing rate 1/r^ that would become dependent on the exponent a. 

4-4 Discussion: Consistency of the calculation of the dephasing rate and the 
transition width 

Having completed the calculation of the dephasing rate and of the localization 
transition width, we have to come back to the assumptions made in course of 
the calculation and check their consistency. Specifically, we have made two im- 
portant assumptions: on the averaging procedure and on the critical character 
of the wave functions involved. 



4-4-1 Critical character of wave functions 

When we calculated the dephasing rate as the imaginary part of the self energy 
ImE(£,e), we assumed that the zero energy (i.e. the position of the chemical 
potential /x) is exactly at the critical point. On the other hand, the obtained 
transition width scales as oc T K with k < 1, i.e. it is much larger than T. 
Therefore, it is important to check that finite deviation of fj, from criticality, 
|/i — /i c | oc T K does not invalidate the calculation. This detuning from criti- 
cality will produce a characteristic scale (localization or correlation length), 
£ M oc \fi — fi c \~" that determines the range of critical behavior of correlation 
functions. Inserting here \fj, — fi c \ oc T K , we get £ M oc T~ l l ZT , which is just the 
condition £ M = that we used to determine the transition width. Therefore, 
for characteristic detuning from critical energy ~ \fi — fi c \ the critical corre- 
lations extend up to oc T^ 1 ^ ZT . On the other hand, the range of spatial 
integration in Eq. (18) was set by the thermal length Lt oc T~ l l d . We thus 
have to compare and Lt- Since Zt < d, we have ^> Lt, which justifies 
the calculation. 
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4-4-2 Averaging of dephasing rate 



While evaluating r^, we have performed averaging over the disorder realiza- 
tions. A natural question to ask is whether this is a valid procedure. Indeed, 
we know that in the localized phase such an averaging breaks down with low- 
ering temperature [36] since the level spacing of those states with which the 
given state is coupled becomes larger than the averaged dephasing rate. As 
a consequence the Golden-rule calculation of the dephasing rate that treats 
these states essentially as a continuum breaks down. However, in the critical 
regime we are considering the situation is essentially different. As we discussed 
in Sec. 4.4.1, the localization length £ M , which is equal to the dephasing L^, 
increases (with lowering temperature) faster than T~ l / d . As a result, the num- 
ber of states contributing to the dephasing rate of a given state — i.e. of states 
located in the energy interval of width T and in the spatial volume (area for 
d = 2) of extension £ M = oc T~ X I ZT increases as a power law with decreasing 
T. In this situation, the dephasing rate is a self- averaging quantity. 

4-4-3 Other possible contributions to dephasing 

Strictly speaking, what we have calculated in this paper is the lowest-order 
(golden rule) dephasing rate governed by the generic range of frequencies (f2 ~ 
e' ~ T). Thus, putting it rigoristically, we have calculated the upper bound on 
Ltf,, and, consequently, the lower bound on zt and upper bound on k. Indeed, 
formally one cannot exclude the possibility that a contribution of a different 
scaling domain of frequency and spatial variables to the lowest-order dephasing 
rate or contribution of higher order in U is larger. Our preliminary analysis 
indicate that this does not happen (apart from possible logarithmic corrections 
to scaling) for critical points with relatively small anomalous exponents, like 
Anderson transition in 2 + e dimensions and quantum Hall transition which 
are in the focus of this paper. It is a challenging task to understand whether 
such exotic contributions may dominate the dephasing rate at transitions with 
strong fluctuations, like Anderson transition in higher dimensionalities. The 
complexity of the problem is related to the fact that the other contributions 
are controlled by wave function correlations characterized by different scaling 
exponents. We postpone this analysis to future work. 



5 Coulomb interaction and experiment 

Most of this paper is devoted to the case of short-range interaction, when the 
interaction irrelevant in the RG sense and the dephasing rate is controlled 
by critical properties at the non- interacting fixed point. In this section we 
briefly discuss the present understanding of the case of long-range Coulomb 
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interaction and summarize the available experimental results. 

5.1 Criticality at Anderson and quantum Hall transitions with Coulomb in- 
teraction 

For long-range (1/r) Coulomb interaction, the dephasing rate is proportional 
to temperature, ~ T. It is easy to see how the short-range interaction re- 
sults develop into this behavior with decreasing A. Indeed, the lowest possible 
"short-range" A is A = d, for which the first line of Eq. (21) should be used, 
yielding l/r^ ~ T. Further decreasing A does not change the scaling of the 
dephasing rate anymore, since l/r^ cannot vanish slower than temperature in 
a system with meaningfully defined fermionic excitations. The result ~ T 
implies that the scaling with temperature is the same as with frequency, which 
is usually the case at "standard" quantum phase transitions. Therefore, in con- 
trast to the case of short-range interaction, for the long-range interaction there 
is no need in distinguishing the dynamical exponents governing the frequency 
and temperature scaling, zt = z. 

On the other hand, the long-range interaction problem is characterized by 
several dynamical exponents controlling the frequency scaling of different ob- 
servables [3,4]. The reason for this complex behavior is the existence of several 
conserved quantities. Let us assume for simplicity that the system is spin- 
polarized (or else, the spin invariance is completely broken, e.g., by spin-orbit 
interaction or by magnetic impurities). Then the conserved quantities are par- 
ticle number and energy. As a consequence, there are two Goldstone modes, 
which are characterized by poles of the diffusion type but with a non-trivial 
scaling, q Zi ~ to. In notations of Ref. [4] z\ corresponds to the energy mode, 
and Z3 to the density mode. 

The exponent zx, which was denoted as ( by Finkelstein [3] and as 2 + 7* in 
Ref. [39], governs the renormalization of frequency in the cr-model action. It is 
this exponent that plays a role of zt = z for the problem we are considering, 
i.e., it controls the scaling of the dephasing length ~ tJ z and therefore 
enters the formula 

k — 1/vz] z = Z\ (65) 

for the exponent k of the transition width. 

We emphasize the distinction between different dynamical exponents, since 
this was not always appreciated by researchers in the field and has led to con- 
fusions and controversies. A number of authors used the exponent 23, which 
is equal to 1 for the case of 1/r interaction, instead of z\ in Eq. (65). This is 
incorrect. To make this point clear, it is constructive to draw an analogy with 
a 2D problem: the critical problem in 2 + e dimensions with small e bears a 
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lot of similarities with a 2D problem with large conductance. The exponent 
z 3 = 1 exists also in 2D: it controls the plasmon pole in the (reducible) density 
response function. It is well known, however, that the plasmon pole does not 
affect the conductivity; in particular, the dephasing rate is recalculated in de- 
phasing length according to a diffusion formula oc (Dt^) 1 ^ 2 , corresponding 
to z — 2. Furthermore, the localization effects are controlled by cooperons 
and by "delayed diffusons" [37,38] that are given by ladder diagramms with- 
out interaction vertex corrections. These are just modes that acquire the z\ 
dynamical scaling at criticality. 

In 2 + e dimensions for the considered symmetry class (time reversal and spin 
symmetries are broken; denoted as "MI(LR)" in Ref. [4] ), the /3-function, the 
critical point, and the exponents v and z are known up to the two-loop order 
[39]: 



P(t) = et - 2t 2 - 4At 3 ; A ~ 1.64 ; (66) 

U= e --p + 0{J) (67) 

^=1-^ + 0(6); (68) 
e 

- 2 +HMi-3)' 2+0 <' 3 »- (69) 

Substituting (68) and (69) into (65), we get 



For the quantum Hall transition with Coulomb interaction not much can be 
said concerning the values of the exponents on the theoretical level. Since the 
transition happens in the strong-coupling regime of the cr-model, no reliable 
analytical predictions can be made. The numerical analysis is also very diffi- 
cult: (i) exact diagonalization can be performed for small systems only, which 
is not suffucient for determining the critical behavior, and (ii) no approxi- 
mate method that would allow to get the interacting critical exponents has 
been developed. On general grounds, and using the analogy with the Anderson 
transition in 2 + e dimensions, one can say that there is no reasons to expect 
that the exponents v and z (as well as any other exponents) would take the 
same values as for the non-interacting system. Also, there is no reasons for 
the exponent z = z\ to take any "simple" value (like, e.g., 1 or 2). 
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5.2 Experiments on localization transition in interacting systems 

5.2.1 Anderson transition in 3D 

The 3D localization transition was extensively studied on doped semiconduc- 
tor systems, such as Si:P, Si:B, Si:As, Ge:Sb. In most of the works, samples 
with a substantial degree of compensation [i.e. acceptors in addition to donors, 
e.g. Si:(P,B)] were used, which allows one to vary the amount of disorder and 
the electron concentration independently. On these samples, values of the con- 
ductivity exponent s in the vicinity of s ~ 1 were reported [40] with scattering 
of values and the uncertainties of the order of 10%. A similar result was ob- 
tained for an amorphous material Nb x Sii_a;. We recall that s is expected to 
be equal in 3D to the localization length exponent v according to the scaling 
relation s = v[d — 2). 

On the other hand, the early study of the transition in undoped Si:P [41] 
gave an essentially different result, s ~ 0.5, which is also in conflict with the 
Harris inequality v > 2/d. This discrepancy was resolved in [42] where it was 
found that the actual critical region in an uncompensated Si:P is rather narrow 
and that the scaling analysis restricted to this range yields s ~ 1.3. A more 
recent study along these lines [16] yielded s = 1.0 ±0.1, in agreement with the 
values obtained for samples with compensation. It thus appears that for the 
orthogonal symmetry class (preserved time-reversal and spin invariances) the 
experiments have converged to the value s = 1.0 ± 0.1. (The only exception 
is a recent experiment on uncompensated Si:B [17] where a larger value was 
found, s ~ 1.6. A possible explanation is that the temperatures reached in 
this work were not sufficiently low. Another possibility is that Si:B belongs to 
a different universality class, in view of stronger spin-orbit scattering.) Results 
on the dynamical scaling are much scarcer: it was found to be z = 2.94 ± 0.3 
in Ref. [16] and z « 2 in Ref. [17]. 

The fact that the experimental value s ~ 1 differs from the result s = v ~ 
1.57 ± 0.02 [43] for non-interacting systems of the orthogonal symmetry class 
is in line with the general expectation that the Coulomb interaction affects 
the critical exponents. 

5.2.2 Quantum Hall transition 

Experiments on the integer quantum Hall plateau transition determine the 
width of the critical region (peak in a xx and plateau transition in o~ xy ), which 
scales with the temperature T as ABocT K . The early measurements of the 
exponent k performed on InGaAs/InP samples yielded k=0.42±0.04 [14]. A 
number of works discussed the effect of macroscopic inhomogeneities [46,47,15] 
that complicate observation of the true IQH critical behavior. The final con- 
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elusion is [15] that for short-range disorder, when the true IQH criticality can 
be achieved, k=0.42±0.01, in agreement with the result of Ref. [14]. A very 
recent work [48] where the quantum Hall transition in AlGaAs/AlGaAs sam- 
ples was analyzed down to very low temperatures (lmK) confirmed this result 
for k. 



Experimental determination of the exponents v and zj- is a highly complicated 
problem. The values reported in the literature are v ~ 2.3 [44,45] and z? — 1 
[48], with the scattering of data of the order of 10%. Our feeling, however, is 
that these data might be essentially affected by systematic errors. Specifically, 
the works [44,45] where v was measured reported the values of k in the range 
from 0.6 to 0.8 (i.e. much larger than the true exponent 0.42). It is understood 
now that such increased values of k correspond to situations where macroscopic 
inhomogeneities do not allow one to observe the true quantum Hall criticality. 
In such a situation the observed v may also differ from the actual quantum 
Hall value. Further, the most direct determination of z [48] was based on the 
analysis of the data for different system sizes, which is a nice way to find 
the dephasing length. However, there may be a problem [49] related to the 
fact that disorder strength was apparently correlated with the sample width 
(see Fig. 3a of Ref. [48]). As a result, the natural temperature scale (mean 
free time) for different samples is different. This is not a small effect: as is 
seen from Fig. 3a of Ref. [48] changing the sample width by factor of 5 not 
only changes the saturation temperature by factor of ~ 5 but simultaneously 
changes the characteristic temperature scale by factor of ~ 2.5. In our view, 
this might considerably affect the determination of the dynamical exponent 
zt- It seems that more experimental work may be needed to overcome these 
difficulties related to systematic errors in evaluation of zt and v. 



Summarizing the experimental findings, the most updated values of the expo- 
nents are v = 2.3±0.1, «=0.42±0.01, and z T = z ~ 1.0 ± 0.1, although the 
error in determination of v and Zt may be considerably underestimated due 
to systematic errors. Let us remind the reader that the theoretical results for 
the case of short-range interaction are as follows: the value of v ranges from 
2.35 to 2.59, zt — 1.23, and k = 1/vzt is in the range from 0.314 to 0.346. 
It appears that the difference in values of the exponents between the cases of 
short-range (theory) and long-range (experiment) interactions is not so large: 
< 10% for v, < 20% for z T , and < 30% for k. (Again, for v and zt might be 
large due to systematic errors.) Nevertheless, the difference demonstrates that 
the current experiments on criticality at quantum Hall transitions cannot be 
explained in terms of the non-interacting fixed point. An experimental real- 
ization of the short-range interaction universality class remains a challenging 
issue for future research. 
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6 Conclusions 



To summarize, we have studied the scaling properties of dephasing rate at crit- 
ical point of the localization transition. We have considered the case of a short- 
range interaction in systems with no spin degeneracy (or broken spin-rotation 
symmetry). In this situation the interaction is found to be RG- irrelevant, and 
the critical properties can be studied by performing the scaling analysis near 
the non- interacting fixing point. 

More specifically, we considered problems with broken time-reversal invari- 
ance: the quantum Hall transition and the Anderson transition in 2 + e dimen- 
sions. Our work combined analytical and numerical analysis. In the analytical 
part, we used the framework of the non-linear a model. We identified oper- 
ators controlling the scaling of the correlation function that determines the 
dephasing rate. Further, we performed their RG analysis in 2 + e dimensions. 
This allowed us to find the analytical results for the exponents p, zt, and k 
governing the temperature scaling of the dephasing rate, dephasing length, 
and the transition width. 

The numerical analysis was used to obtain critical exponent at the quantum 
Hall transition. Our results for the exponents largely agree with those obtained 
in Refs. [18,19]. However, our system sizes are much largely than those studied 
in Ref. [18] that was crucial for getting a window of distances where corrections 
to power-law scaling are small. 

Experimental results on localization transition with short-range interaction 
are extremely desirable. In the case of quantum Hall transition one could 
imagine screening of the Coulomb interaction by an external gate. It would be 
extremely interesting to observe the change of the exponents compared to the 
case of long-range interaction. This is a challenging task, especially since the 
difference between the exponents appears to be not so large (see Sec. 5.2.2). 
For 3D Anderson transition, its experimental realization and investigation in 
cold-atoms systems would be of great importance. 

On the theoretical side, the wave function correlation exponents \xi and a need 
to be evaluated for 3D Anderson transition (with and without time-reversal 
invariance), in order to predict the behavior of the dephasing rate and thus the 
exponents p, z?, and k. More generally, investigation of the statistics of wave 
functions at criticality beyond the leading multifractal behavior represents an 
important research field. It would be interesting to understand the evolution of 
this statistics (which includes the statistics of Hartree-Fock matrix elements) 
from the regime of weak to strong multifractality, as it has been done for the 
spectrum of leading multifractal exponents A q [1]. 
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A Correlation function K\. 

In this Appendix we demonstrate that the function /Q corresponds to the 
eigenoperator P^i. The function /Ci can be expressed in terms of the exact 
single-particle Green's functions G^^r, r'): 



K x = — ( Im Gf+Jn, n) Im G|(r 2 , r 2 ) - Im G§ + Jr 2 , n) Im Gf (n, r 2 ) 



Kl = To — ^2 (^AQaa{ri)tTAQ bb (r 2 ) + trAQ ab ( ri )AQ ba (r 2 )) , (A.2) 



where repica indices a, b are different, a ^ b, which reflects the fact that we 
are dealing with two different eigenstates. Since we are interested in the case 
for which two points are close to each other (say, | r*i — T*2 1 is of order of a few 
lattice constants a), we can consider the arguments of Q matrices as equal. 
Below the arguments are omitted. 

Let us define two operators bilinear in Q: 




(A.l) 



Following standard steps (see e.g. Ref. [25]), we obtain 



A 2 



0± [Q] = tr AQ ab AQ ba ± tr AQ ab tr AQ ba . 



(A.3) 



Clearly, /Ci oc + . In order to obtain the U(n) x U(n) invariant expression let 
us perform the following global rotation in the Q-matrix space, 



Q(r) U^Q(r)U, ug = V*,** , (A.4) 
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\r n 1/ _ 1 

Vl > 1 ~ n(n 2 -l)> V2 ~ n(n 2 -l)' 



t/ _ n 2 — 2 T 7" _ -it _ 

"1,1,1 — 7n77Z~r\u7Z~A \ i "2,1 — „(W-*-^(r,^-A\ i "3 — 



n(n 2 -l)(n 2 -4) ' " i > 1 n(n 2 -l)(n 2 -4) ' v A n(n 2 -l)(n 2 -4) ' 

t /" _ 6-8n 2 +n 4 T/ _ n(4-n 2 ) T/ _ 6+n 2 

"1,1,1,1 — „ 2 (n 2 -l)(n 2 -4)(n 2 -9)' V2 > 1 ' 1 ~~ n 2 (n 2 -l)(n 2 -4)(n 2 -9) ' V2 2 ~ n 2 (n 2 -l)(n 2 -4)(n 2 -9) ' 

t / 2n 2 — 3 t/ 5n 

V3 < 1 ~ n 2 (n 2 -l)(n 2 -4)(n 2 -9)' ^ 4 ~~ n 2 (n 2 -l)(n 2 -4)(n 2 -9) 

Table A.l 

Coefficients V} for averaging over the unitary group U{n). 



which does not change the action S a . Thus, we introduce the averaged, oper- 
ators 



± [Q] = (trA[U- 1 QU] ab A[U- 1 QU] ba ±trA[U- 1 QU} ab trA[U- 1 QU} ba/i 

(A.5) 

where (■ ■ - )u denotes averaging over U(n) x U(n) global rotations. Since the 
action S a is invariant under U{n) x U(n), the cr-model averages of 0±[Q] and 
0±[<2] are equal, 

(6±[Q\) = (0±[Q\), (A.6) 



In order to perform the averaging over U(n) x U(n) rotations, we use the 
following results [50]: 

((^Wp = V 1 8 ab 8 a ^ (A.7) 

{{U'^loPlbiU^fc^ld) U = Vl,l [SabSa^cdS^ + 5 ad 5 afJ ,5 bc 5^} 

+ V 2 [ShcSapSdaS^ + SahSa^cdSfr] , (A.8) 

where V\, Vi t \ and Vi are given in Table A.l. The index p = ± distinguishes 
between the retarded and advanced sectors, i.e. (U + ,U~) is an element of 
U(n) x U(n). The averaging over U + and U~ is carried out independently. 

Performing the averaging, we obtain the following U(n) x U(n) invariant re- 
sults: 

It is not difficult to check that operators 0± are eigenoperators under the 
action of the renormalization group: 0_ = P2 and + = P\ \ . 
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B Correlation function /C2. 



In this Appendix, we express the correlation function /C 2 in terms of basis 
operators This requires, in addition to bilinear operators considered in 

Appendix A, introducing operators of the fourth order in Q. As in the previous 
Appendix, it is convenient to perform global transformation (A. 4) and average 
over U(n) x U(n) rotations. Then }C 2 becomes 

k^^Irm + UQ]) (b.i) 

where 

Rj [Q] = (tr k(U- l QU) ab k(U~ l QU) ba tr k(U~ l QU) cd k(U~ l QU) dc ) U} (B.2) 

R 2 [Q] = (tr A(t/- 1 g[/) Q6 A(f/- 1 gt/) fec A(^ 1 gf/) cd A(f/- 1 gf/) da )j / .. (B.3) 

In order to perform averaging over U (n) x U(n) rotations it is convenient to 
use diagrammatic technique developed in Ref. [51]. The result is as follows 



Ri[Q] = E V 2<2 (TiQ p ) 4 + (W 4 + 2V 2 ,i,i) TtQ 2 p (TtQ p ) 2 
p=± ^ 

+8V 3 ,i Tr Ql Tr Q p + (2V 2<2 + Tr Q 2 p Tr Q\ 

+(4y 2l i,i + 2V 4 ) Tr Ql - W 1 V 2>1 (Tr Q p ) 2 Tr A p 
-WxVx,!,! Tr Q 2 p Tr A p - 8X4^2,1 Tr(A^) 
-8V 1 V 3 Tr Q p Tr A P Q P + 2(V? A + V 2 2 )(Tr A p f 
+W 1A V 2 Tr A 2 p + V? tl Tr Q 2 p Tr Q 2 _ p 

+2V 2 V ltl TrQl(TrQ_ P Y + V 2 2 (TrQ p ) 2 (TrQ_ p ) 2 }, (B.4) 



MQ] = E ^(Tr g p ) 4 + (4V3,i + 2^2) TrgJ(Tr g p ) 2 
p=± ^ 

+4(F 4 + F 2 ,i,i) Tr Ql Tr g p + (V4 + 2V 2 ,i,i) Tr gj Tr Q 2 p 
+ (V 2i2 + + 4V 3 ,i) Tr gj - 4W 3 (Tr g p ) 2 Tr A p 

-41/^2,1 Tr gjTr A p - 4Vi(F 3 + Tr(A p gJ) 
-8ViV 2l i Tr g p Tr A p g p + 2V 2 2 Tr g p Tr g_ p Tr A p 
+4 VWi,i Tr Q p Tr 5 P + 2V; 2 ! Tr B P Q P + 2F 1)1 F 2 (Tr A p ) 2 

+ (y i 2 1 + K 2 )TrA 2 ], (B.5) 
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where the coefficents Vj lt ,,,j m are introduced in analogy with Eqs. (A. 7), (A. 8). 
Specifically, Vj u ,„ t j m arises when one averages a product of j = j\ + . . . + j m 
matrix elements of U and j matrix elements of U~ x as a coefficients in front 
of terms corresponding to m "cycles" of the sizes ji, . . .j m . Further, we have 
introduced the notations 



tpA 1+pA ( . 

Qp = — 2 — ^ — s — ' ^ ' 



A^-^Q^Q 1 -^, (B.7) 
Bp= l+M e l^A e l^A e l^A (B8) 

Using the expressions for the coefficents from Table A.l, we can express the 
operators R\ and R 2 via operators (29)-(35). The result is given in Eq. (28). 
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